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Abstract 


We examine the decays of excited heavy mesons, including the leading power 
corrections to the heavy quark limit. We find a new and natural explanation 
for the large deviation of the width of the D i(2420) from the heavy quark 
symmetry prediction. Our formalism leads to detailed predictions for the 
properties of the excited bottom mesons, some of which recently have been 
observed. Finally, we present a detailed analysis of the effect of power cor¬ 
rections and finite meson widths on the angular distributions which may be 
measured in heavy meson decays. 
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I. INTRODUCTION 


The excitation spectrum of charmed and bottom mesons has received considerable recent 
theoretical and experimental attention. The discovery of the Bi and the B 2 mesons [|lj, 
and the measurement of their masses and widths, complements the improving data being 
acquired on their charmed cousins, the D i and D 2 p[. At the same time, the theoretical 
understanding of the production and decay of these mesons has profited from the application 
of the Heavy Quark Effective Theory (HQET) and the ideas behind it, in particular the 
enlarged spin-flavor symmetry of QCD which obtains in the limit m c ,mft —> oo 0. 

We will begin this paper with a review of the experimental situation, and of the im¬ 
plications of heavy quark symmetry for excited heavy mesons. We will then investigate 
systematically the effect of the leading corrections to the heavy quark limit. We focus on 
corrections which violate the heavy spin symmetry, since it is predictions which follow from 
this symmetry which are tested by current data. We will find a new and natural explanation 
for the anomalously large width of the D\. Our formalism leads to detailed predictions for 
the properties of the Bi and B 2 , as well as for their strange counterparts. We close by 
examining various angular distributions in the strong decays of heavy mesons, since with 
sufficiently precise data they will eventually provide detailed tests of our proposal. 


II. THE HEAVY-LIGHT MESON SPECTRUM 

The excitation spectrum of heavy mesons takes a particularly simple form in the limit 
m —> oo, where m = m c or m = mb- In this limit, the spin of the heavy quark decouples, and 
both the spin J of the meson and the angular momentum J \ of the light degrees of freedom 
become good quantum numbers. For each state in which the light degrees of freedom have 
spin-parity J[ , there is a degenerate doublet of meson states with J p = J[ ± The mass 
splitting between these doublets arises only from effects of order Aq CD /m in the heavy quark 
expansion. 

In this language, the ground state heavy mesons M and M* have light degrees of freedom 
with J[ = | , and there are low-lying excited P- wave states with J[ = | + and J[ = | + . 
The current experimental situation is summarized in Table [I]. Here we quote errors only 
for the excited D mesons [|2]||]. The data on the excited B mesons is the result of a fit for 
which no errors are given [1J, and the different charge states of the B*, Bi and B 2 have not 
been resolved. States for which no masses are given in Table [I] have not yet been observed, 
although they are expected to exist. 

Heavy quark symmetry imposes a number of constraints on the strong decays of these 
states. Since such decays are entirely transitions of the light degrees of freedom, the four 
possible decays of the two members of a doublet with given J[ to the two members of 
another doublet with J[' are all essentially a single process. In the strict m —> oo limit, in 
which all such doublets are completely degenerate, this fact leads to the simple prediction 
that the two excited states should have exactly the same width. In fact, it is more accurate 
to use the actual masses of the states in calculating phase space effects, imposing the heavy 
quark symmetry only on the level of the matrix elements. Even though this approach is not 
technically consistent in the sense of the 1/m expansion, it allows us to incorporate certain 
1/m effects which are numerically quite substantial. 
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TABLE I. 


For example, since the decay ~^ Dn can take place only with the n in an orbital 
d- wave, the same also must be true of the transitions —■> D*tt and D\ —> D*tt (the 

matrix element for D\ —> Du vanishes). Incorporating phase space effects, the partial width 
for the transition Dj —> Dfir is proportional to , times a mass-independent 

matrix element. This fact, plus some elementary spin counting, allows one to make the 
predictions || 


no? -> £^n = „, r(c°) 

r(B5° D‘+ir-) ' ’ r(ZJJ°) 


( 2 . 1 ) 


(We choose the modes for which the data are most accurate.) Experimentally, one finds [!|.|4|][] 


r(L>*° -► D+tt~) 
T(Df -»• D*+ 7T-) 


2.2 ±0.9, 


np _?) 

r (Df) 


0.71. 


( 2 . 2 ) 


Clearly, one of these predictions works extremely well, while the other works not at all. Why 
might this be so? 

One common explanation is that the D\ has a small mixing with the D[, which decays 
in an s-wave rather than a d- wave, and hence is expected to be considerably broader ||. 
Such a mixing is allowed when spin symmetry violating 1/m effects are included, as the 
Di and D[ both transform as J p = 1 + under the Lorentz group and differ only in their 


x We extract the ratio T(Zl] ) )/r(Zl 2 0 ) from the data in Table |. We do not assign an error, because 
the correct error depends crucially on the correlations between the measurements of the widths, 
which we do not know. 
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values of J[. Since even a small mixing is important if it is with a much broader state, this 
provides a simple explanation of why a 1/m correction of the natural size might lead to an 
anomalously large correction to the total width of the D\. 

Unfortunately, it is not one which is particularly favored by the present data. One may 
measure the angular distribution of the emitted pion and determine directly whether it is 
in an s-wave or a d- wave. The situation is complicated by the fact that the pion angular 
distribution depends not only the ratio of the s and d -wave partial widths, but also on the 
relative phase 77 between the two matrix elements. Particularly in their data on D® decay 0 , 
CLEO finds that a large s-wave component is compatible only with a restricted region in 
cos 77. By no means is such a scenario ruled out, but it is less than generic. 

In addition, there is no evidence for a significant s-wave component in the decay of 
the D s 1 , which is related by flavor SU(3) symmetry to the D\. One may use heavy quark 
symmetry to predict the d -wave width of the D s i in terms of the width of the D* 2l analogous 


to the second relation of Eq. (O). The D s \ and D* 2 decay via K emission to D and D* 


Given the measured T(.D* 2 ), one predicts 


T(D S 1 ) = 0.3 MeV, 


(2.3) 


far below the CLEO upper limit. Still, there is little room for a large additional s-wave 
component. To see this, one may use flavor SU(3 ) and the upper limit on T(ZUi) to predict 
an upper limit on T(T ) i) | 6 |]. Since || in D\ —> D*n is typically larger than p^ in 
D s 1 —> D*K , the upper limit on T(.Di) will be much more stringent if the decay is assumed 
to be primarily s-wave, as opposed to d- wave. One then obtains independent upper limits 
on the s-wave and d -wave components of T(Z^i): 

T s (Di) < 3MeV, (2.4) 

Td(i4 1 ) < 105MeV. 

The correlated limits are somewhat stronger (see Ref. [(D). However, we see immediately 
that under these assumptions it is impossible to accommodate a large s-wave component in 
D\ decay. 

Of course, these assumptions might not be very good. Flavor SU(3) could fail badly here. 
One test of SU(3) is to use the width of the D 2 to predict the width of the D* 2 , assuming d- 
wave decay and including the correct phase space for the K. One finds T(Z1* 2 ) = (9±3) MeV, 
in reasonable agreement with experiment. (Perhaps it would be more correct to include the 
577(3) violating factor /^//| in this prediction; however, doing so makes the agreement with 
experiment worse.) However, it is possible that the D s i — D' sl mixing is very different from 
the D\ — D\ mixing, since the angle depends delicately on the interplay of a mixing matrix 
element and a mass splitting, both of which receive 5f/(3)-violating corrections. Still, we 
are not encouraged that this explanation for the anomalously large Di width is the correct 
one. 

As an alternative, it has been suggested 0 that the Di width receives a large contribu¬ 
tion from the emission of two pions which resonate through a p meson, with no analogous 
enhancement in D 2 decay. Two pion decays in which one of the pions resonates with a broad 
D\ or Dq have also been considered j5|, but are not thought to contribute significantly to 
the total width. In the next section, we will introduce a simpler explanation for the D 1 
width, which arises naturally at higher order in the heavy quark expansion. 
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III. THE CHIRAL LAGRANGIAN 


The strong decays of excited mesons involve the emission of soft pions and kaons, and 
hence it is useful to analyze these interactions with the help of chiral perturbation theory. 
To be concrete, we will specify to the charm system; the generalization to bottom is at all 
points straightforward. The chiral lagrangian appropriate to the analysis of ground state 
and excited heavy mesons has been derived elsewhere [ 0 ; here we will simply recall the basic 
points. 

The heavy mesons are represented by matrix superfields which carry a representation not 
only of the Lorentz group but also of the 577(2) heavy quark spin symmetry. In addition, 
they transform as 3’s under flavor S77(3), since by convention our heavy mesons contain 
a single heavy quark (rather than an antiquark). For the ground state and lowest excited 
states discussed in the previous section, these superfields take the form [ITU 


H a = 


S n = 


± a 


(i + i>) 

2 V 2 

(1 + 1O 

2 V 2 
(i + fl 
2V2 


- Dorr 


DTa 7m7 6 - D< 


0 a 


(3,1) 


- d 


la 


7 5 ( 5 ; 




where v M is the four-velocity of the heavy meson, and a is the flavor index. These fields 
are normalized nonrelativistically. Under a heavy quark spin rotation Sq, the superfields 
transform as 77 —> SqH , etc., while under a Lorentz transformation Sl they transform as 
77 -> 57.7751, etc.. 

The heavy mesons interact with the octet of pseudo-Goldstone bosons, which are treated 
with the usual nonlinear formalism. The lagrangian is written in terms of an exponentiated 
matrix of boson fields, £ = exp(LM// 7 ), where 


/ , /I -0 


M = 


7T U + 
7 T~ 

K~ 


I 7 


7 r T 

\ 7r ° + 
77° 


K+ \ 
K° 

Vj 


(3.2) 


and f n ~ 135 MeV. Under SU(3)l x SU(3)r, the field £ transforms as £ —» L£ W = 7/£7?f, 
where U is a matrix which depends on Ai. For generators within the diagonal subgroup, 
U = L = R. 

The pseudo-Goldstone bosons couple to the heavy superfields through the covariant 
derivative D^ b = 5 a bd^+\(7d^i-\-id^7)ab and the axial vector held A% b = f (£ t 5 /i £—£cU£la&- 
Under SU(3) L x SU(3)r, the helds transform as A -> UAU\ 77 -»■ 7/77, 77^77 -> 7/77^77, 
and so on. The kinetic part of the chiral lagrangian then given by 


£ kin = -/ w 2 Tr 


d^d u T) 


Tr 


771/7 


-Tr 


Hiv ■ DH 


+ Tr 


S(iv ■ 77 - A s )S 


T (iv ■ 77 - A t )T^ 


(3.3) 


where S = £ 2 . The trace is taken with respect to spinor and flavor indices, which we 
suppress. The excitation energies A 5 and At are defined with respect to the spin-averaged 
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masses of the excited doublets. The heavy fields obey the equations of motion iv ■ DH = 0, 
iv ■ DS = A s S and iv ■ DT M = A 

The interactions between the various fields are constrained by heavy spin symmetry 
and chiral symmetry, ^ as well as by Lorentz invariance (including parity and time reversal) 
and velocity reparameterization invariance 1 ( 113 , 0 . We are interested in matrix elements 


between the ground state doublet and each of the excited states, for which we introduce the 
interaction terms 


C a = f Tr iHS^A^+Kc., 


c -‘ = r Tr 


HT> i 'fy i {iD li A v + iD v A li ) 


+ h.c. 


(3.4) 


These are the leading terms in the chiral expansion, which is an expansion in derivatives and 
fields over the chiral symmetry breaking scale A x ~ 1 GeV. Note that Cd, which mediates 
the d- wave decays of the D\ and D appears at higher order than C s , which mediates the 
s-wave decays of the Do and D[. Velocity reparameterization invariance plays a particularly 
important role in constraining C ( ] [|6j]. 

One may include bottom mesons by introducing new superfields H b , S b and Tjf to repre¬ 
sent the ground state and excited B mesons. These fields appear in analogues of the terms 
Akin, C s and Cd, with the same parameters A s, A t, f and h as in the case of charm. This 
doubling of terms persists at higher orders, again with identical coefficients, except that 
explicit factors of 1/m should be taken to be l/m c or l/m b as appropriate. 

The chiral lagrangian, as developed so far, is sufficient to make the heavy spin symmetry 
predictions discussed in Section II. However, as we have seen, certain of these predictions 
work better than others; in particular, the ratio of partial widths —■► D*n) /T(Z ?2 — 5 ► 

D7 r) is well predicted by the heavy spin symmetry, while the ratio of full widths T(£)i)/T(£)|) 
is not. Evidently, corrections to the heavy quark limit can be quite large. Why is this so 
in some cases, but not in others? To gain insight into this question, it is necessary to go to 
subleading order in the 1/m expansion. 

The most general extension of the chiral lagrangian to order 1/m is extremely unwieldy. 
However, we are interested only in those 1/m corrections which break the heavy spin sym¬ 
metry. We will find that when we neglect systematically spin symmetry conserving terms, 
the resulting theory is sufficiently constrained to yield interesting information. 

The 1/m corrections to the chiral lagrangian arise from 1/m corrections to the heavy 
quark lagrangian, which is derived from the full QCD lagrangian in the m —> 00 limit. This 

m 


expansion takes the form 


C 


1 


HQET 


= hiv ■ Dh + -— h{iD) 2 h + -— ha^ v (\gG + 


2m 


2m 


(3.5) 


where h is the HQET heavy quark held. The effect of the subleading terms h(iD) 2 h and 
gha^G^h on the chiral expansion may be treated in the same manner as other symmetry 


2 Since we are only interested in transitions involving a single light meson, flavor SU( 3), under 
which the pions transform linearly, would be sufficient for our analysis. We employ the full formal¬ 
ism of nonlinear representations simply as a convenience. 
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breaking perturbations to the fundamental theory such as finite light quark masses. Namely, 
we introduce a “spurion” held which carries the same representation of the symmetry group 
as does the perturbation in the fundamental theory, and then include this spurion in the 
chiral lagrangian in the most general symmetry-conserving way. When the spurion is set to 
the constant value which it has in QCD, the symmetry breaking is transmitted to the effective 
theory. In the case of finite light quark masses, for example, the symmetry breaking term in 
QCD is qM q q , where M q = diag (m u , rrid, m s ). Introducing a spurion M q which transforms as 
M q —> LM q R f under chiral 577(3), we then include terms in the ordinary chiral lagrangian 
such as /j Tr [M q Y> + M q TX], 

In the present case, only the second of the two correction terms in £ H qet violates the 
heavy spin symmetry. We include its effect in the chiral lagrangian by introducing a spurion 
which transforms as —> Sq^ v Sq under a heavy quark spin rotation Sq. This 
spurion is introduced in the most general manner consistent with heavy quark symmetry, 
and is then set to the constant to yield the leading spin symmetry violating 

corrections to the chiral lagrangian. We will restrict ourselves to terms in which appears 
exactly once. 

The simplest spin symmetry violating effect is to break the degeneracy of the heavy 
meson doublets. This occurs through the terms 

Xh Tr - A 5 Tr \s$T'Sa lu/ ] - A T Tr [r"$f . (3.6) 

The dimensionful coefficients are fixed once the masses of the mesons are known. For the 
ground state D and D *, for example, we find 

X H = l \Mq* - Ml] = (260 MeV) 2 . (3.7) 

8 L 

This value is entirely consistent with what one would obtain, instead, with the B and B* 
mesons. For the Di and D%, we find 

A T = A \m 2 d * - Ml 1 = (190 MeV) 2 . (3.8) 

16 L 2 1 J 

Note that \/\h and a/At are of order hundreds of MeV, the scale of the strong interactions. 
We are interested in the spin symmetry violating corrections to transitions in the class 
—> Hn, which will arise from terms analogous to C f j, but with one occurrence of 
The spin symmetry, along with the symmetries which constrained Cd, requires that any such 
term be of the generic form 

—— Tr \H^T a C, vaPK7 5 [iD^A K + iD K A + h.c., (3.9) 

where C M ^ a/ 3 K is an arbitrary product of Dirac matrices and may depend on the four-velocity 
v x . This would seem to allow for a lot of freedom, but it turns out that there is only a single 
term which respects both parity and time reversal invariance: 

Cdi = ^- Tr (iD a A K + iD K A a )\ + h.c.. (3.10) 
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We expect the new coefficient hi, which has mass dimension one, to be of order hundreds of 
MeV. 

Finally, there may be additional correction terms which come about by the application 
of velocity reparameterization invariance (VRI) |ll| to the leading interaction term Cd- This 
is a “symmetry’’^ which, as does a gauge symmetry, arises because of a redundancy in the 
variables appearing in the lagrangian. The four-velocity v x which describes the heavy meson 
held is arbitrary up to terms of order 1/m, and the lagrangian must be constructed so as to 
be invariant under reparameterizations of the form v x —> v x + e x /m, where the components 
of e A are of order Aqcd- The heavy meson fields also transform nontrivially under velocity 
reparameterization. VRI is a symmetry which constrains the new terms which may appear 
at higher order in the 1/m expansion in terms of those which are already there. The terms 
which VRI generates at order l/m n+1 may be found by making the replacements §0 


w A -» v x H- iD x , 

m 

H ^ H + ^ Y ' iD - H] ’ 
s ^ S + ,ID U S} , 

-|—— [Y, iD u T^] - —v fl iD u T p 
2 m m 


(3.11) 


in the lagrangian at order 1/m. By the same token, all terms at order l/m n+1 with deriva¬ 
tives acting on heavy meson fields must be consistent with such replacements at one order 
lower. 

New interaction terms of the same order as Cdi will be generated when we make the 
replacements (|3.11|) in the leading term Cd■ However, massaging the new terms with inte¬ 
gration by parts and application of the equations of motion, we fold that they all may be 
written in the form Tr HT^f^d, A) , for some Dirac matrix valued function f(d, A). Hence 
they do not break the heavy spin symmetry and we may ignore them for our analysis. 

The other potentially important effect of 1/m corrections on the decay of the Di is a 
possible mixing between this state and the D[. The D\ transforms the same as the Di 
under the Lorentz group, but differently under the heavy spin symmetry, so spin symmetry 
violating effects can mix the two states. Such a mixing can have a dramatic effect on the 
width of the D\, which now may decay via s-wave pion emission. 

The following parity and time reversal invariant term will induce a nonzero matrix ele¬ 
ment for the mixing of the Di and D\: 


Cmix — .91 Tl 




+ h.c. 


(3.12) 


Unfortunately, the magnitude of the mixing matrix element does not by itself determine the 
mixing angle -0. To find 0, we also need the lowest order splitting A j. — A $ between the 
masses of the two states. Then the mixing angle is given by 


3 We use quotes because this is not a physical symmetry of nature, i.e. there is no associated 
conserved charge. 


















(3.13) 


tan 'ip 


6 , 


where 5 = |(A t — As) and 8 g = —\j\g\jm. Since the D[ has not yet been observed, [] we 
do not know As; hence it is more convenient to treat ^ itself as a free parameter, rather 
than write it as in Eq. (j'i . 1 3|) . 


IV. EXPERIMENTAL IMPLICATIONS 


A. The D\ and D 2 widths 


The new interaction term affects the decays of Di and D 2 in ways that do not 
necessarily respect the heavy spin symmetry predictions 0 . It is straightforward to 
compute the single pion partial widths of these excited states in terms of the coupling 
constants h and h\. Lorentz invariance requires that the decays of the D 2 still involve the 
emission of a d -wave pion, and we find 


T(D* 2 ° -> Dtt) 
T{D° 2 * -> D*tt) 


1 M d 4\p n \ 5 
IOtt M d * A2/2 

3 M D * 41 | 5 

20 vr M d - 



m c 



m c 


(4.1) 


Here and in Eq. (|4.2|) below we include emission of both charged and neutral pions, neglecting 
the small phase space differences between the two channels. Note that the heavy quark 
symmetry prediction for the ratio of the D 2 partial widths is unaffected by the correction 
C ( j,\ ■ This is good, because we have seen that the lowest order prediction works quite well 
already. Truncating Eq. (|4.1|) at lowest order, inserting the experimental width of the D 2 , 
and taking A x = 1 GeV, we obtain the estimate h ~ 0.3. 

The effect of C ( i\ on D i decay is more complicated, because the new term can mediate 
both d -wave and s-wave decays. The decay width is given by 


r(£>i 


D*7t) = — 


1 M D * 4|p„ 


4vr M Dx A 


h + 


5hi 
3 m r 


+ 


8 h\ 

9 rrii 


(4.2) 


The Erst term corresponds to a d -wave pion and the second to an s-wave pion. Note that 
here the s-wave width is also suppressed by || 5 , and is in no sense intrinsically larger than 
the d -wave width. It is consistent to neglect it at this order, since we have not included 
1/m 2 interaction terms in the lagrangian. The corrected ratio of the D 2 width to the Di 
width is then given by 


4 It is possible that the Dj( 2440)^, reported by the TPC Collaboration [13] in the D *°7r :t channel 
with a width of 40MeV, is the D'±. However, this observation has not been confirmed. 
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r (d°i) 
r (Df) 


0.30- 


(4.3) 


1 + 


16 hi 

3 m c h 


+ o 



where we assume that the widths are saturated by the one pion decays. We see that in 
this prediction, the correction is enhanced by a large numerical prefactor. If we generalize 
Eq. (|4.3|) to include a mixture of d -wave and s -wave pion emission, we find 


rpi) 

r(£>|) 


0.30 


cos 2 0 



16/lj \ 

3 771 c h ) 


+ sin 2 0 


3 P 


El 


4(VAJ 2 IPt, 


(4.4) 


With A x = 1 GeV, h = 0.3, and m c = 1.5 GeV, this reduces to 


rpi) 

T(D* 2 ) 


0.30 



h! \ 
85 MeV ) 


cos 2 0 + 77 f 2 sin 2 0 


e ^' 0.71. 


(4.5) 


We may saturate the experimental width of the D\ with d- wave decay by taking h\ = 
115 MeV, which is quite reasonably small in view of the sizes of the similar corrections A h 
and \t■ Alternatively, if hi = 0, then (77f 2 — l)sin 2 0 = 1.37, or / sin0 « 0.13 if / is of 
order one. 


B. Predictions for excited B mesons 

The heavy spin symmetry may also be combined with the heavy flavor symmetry to 
predict the masses and widths of the B mesons in terms of those in the charm sector. The 
splittings between excited doublets and the ground state should be independent of heavy 
quark mass, while spin symmetry violating intradoublet splittings scale like 1/m. If we 
define the spin-averaged masses and mass splittings 

Mb = | Mb* + \Mb , 

Mb** = | Mb* + § Mb 1 , ( 4 - 6 ) 

A Mb** = Mb* — Mb 1 , 

and analogously for charm, the heavy quark symmetries predict 

Mb** — Mb = Mb** — Mb , 

772 

AM B ** = —■ A Mb** ■ (4-7) 

m h 

With m c /mfe = 1/3, and averaging over the charged and neutral charmed mesons, these 
relations yield 


M Bl = 5780 MeV, 
Mb* = 5794 MeV, 
M Bal = 5886 MeV, 
M b * = 5899 MeV. 

■ l s 2 


(4.8) 
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Compared with the data in Table 1|, we see that the measured masses are somewhat lower 
than expected, especially for the nonstrange mesons. The leading corrections to the predic¬ 
tions for Mb** and Mb** are of order 

*~ A «-(i- 2 ^)~ 4 ° MeV ’ < 4 ' 9) 

where we have estimated a QCD scale Aqcd ~ 400 MeV, so the accuracy with which the 
predictions ([4.8Q work is more or less what one would expect. 

ft is somewhat more delicate to make predictions for the widths of the excited B mesons, 
since these depend on the available phase space, hence on the values of the heavy meson 
masses. For the B\ and B s i, they also depend on what one assumes about = | + and 4 + 
mixing in the bottom sector. Let us introduce a notation for the pion momenta which arise 
in these decays: 


B* 2 - 

-> Bn : 

P2 B , 

b;- 

* B*n : 

P-2B* 

B x - 

♦ B*n : 

Pi B* 


and similarly for charm. Then, assuming the dominance of the one pion decay channel, the 
width of the B% is related to that of the D\ via 

T(B*) _ Mp* j 0 A\p 2 b\ 5 Mb + 0-6|p 2 _B* \ 5 M b * \ 
r(I?2) Mb* \0A\p2d\ 5 M d + 0.6\p2D*\'Mp* J 

from which we find the prediction 


r(B%) = (16 ±6) MeV, 


(4.11) 


with the masses given in Table ||. This width is somewhat low, but it is extremely sensitive 
to the mass of the B\ and grows rapidly with Mb* ■ Perhaps here we have a hint that when 
the masses of the Bi and Bi, are better measured, they will be closer to the values predicted 
by heavy quark symmetry. 

We may also generalize Eq. ([4.4| ) to predict the ratio of the widths of the B i and the iLJ, 
assuming once again that one pion decays dominate. Leaving the dependence on the meson 
masses explicit, we find 


£Cgi) 

r(BJ) 


M b * 


Pi B* 

5 M b * 

M Bl 

0.4 \p 2 B\ 5 Mb + 0.6 P 2 b*\ 5 M b * _ 


X 


COS 2 Ipb 



16/ii \ 

3 rribh ) 


+ sin 2 0 6 


3 f 2 E\ b , 1 

4(VAx) 2 |Pis*l 4 J ’ 


(4.12) 


Since the mixing is generated by the spin symmetry violating operator £ rnix ( p,12|) , it should 
scale inversely with the heavy quark mass. Hence for small mixing angles, we might expect 
-0b ~ (m c /nib )^. However, as we see in Eq. (|3.13|) , "06 also depends delicately on the mass 
splitting between the Bi and the B[, so we should put no particular trust in this estimate of 
ijjb- Instead, we will make specific predictions of T(Bi) only in the two limits of pure d- wave 
and pure s-wave decays of both the Di and the Bi, using the masses in Table 
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Y(Bi)/Y(B*) = 0.9 
y(Bi)/y(b* 2 ) = 1.4 


pure d- wave, 
pure s-wave. 


(4.13) 


For the s-wave case, we take ^ = 0 = 7 t/2 and choose / to give the correct Di width. 
While this extreme limit is not favored by the data on D\ decay, it yields a useful upper 
bound on r(T>i). We see that d -wave dominance is somewhat favored by the current data, 
from which one finds Y(Bi)/Y(B 2 ) ss 0.8. 

Finally, we make predictions for the one kaon widths of the excited strange B mesons. 
Applying the analogue of Eq. ( |4.10|) , we find 


Y(B* s2 ) = (7 ± 3) MeV 

(4.14) 

and from the analogue of Eq. (|4. 1 2|) , 


r(s,i)/r(s; 2 ) = o.4, 

(4.15) 


where the B s i is assumed to decay via the emission of a K in a d- wave. These predictions 
will be tested as the data on the B* 2 and B s i improve. 


C. Angular distributions 


While the explanation of the width of the D\ which we have presented is certainly 
consistent with the data, we would like to be able to test it in somewhat more detail. We 
may do so by considering the angular distributions of the pions emitted in its decay. These 
distributions will depend on the chiral lagrangian parameters h, hi and /, and on the mixing 
angle 0 . Finally, with our phase conventions the coefficient of C s may be complex, so we 
will take / —> / exp (A/). 

In addition to the two constraints on the set of parameters {h, hi, /, r], 0 } from the 
experimental widths of the D 2 and the D \, we will now assume that the mixing angle 0 is 
small, ffence we will drop terms which are suppressed by sin 0 , unless they are enhanced 
by a large phase space factor. Note that for values of / of order unity, the constraint ( |4.5|) 
requires that 0 not be too large. An estimate of 0 based on a quark wavefunction model 
gives 0 ~ 9° [T|, which perhaps also supports the use of this approximation. 


1. Two-pion distributions 

In the decays (Di, D 2 ) —> D*tii —> Dttitt 2 ., the angle between the two pions contains 
information about the initial spin state. Let a be the angle between the momenta p^i 
and p,^, as measured in the rest frame of the excited meson. Since the Di and D 2 are 
separated by approximately 40 MeV and have intrinsic widths of 20 — 30 MeV, they overlap 
considerably. The distribution in cos a is a function of where the pion’s energy places it 
in relation to the two resonances. The formalism of Ref. [jnj may be used to extend the 
results of Ref. || to the case where the finite widths Y and Yn- are taken into account. 
We obtain 
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cosa 


(a) 




cosa 

(b) 



(c) 


(d) 


E= 380 MeV E =400MeV E = 420 MeV 

Jt K K 

FIG. 1. The differential distribution in cosa, for the D\ decaying (a) in a pure d- wave; (b) 
in a pure s-wave; (c) in a mixture with h\=0, ip = 9° and / = —0.85; (d) in a mixture with 
hi=0, i\) = 9° and / = 0.85. The solid curve is for E n = 380 MeV, the long dashed curve for 
E n = 400 MeV, and the short dashed curve for E n = 420 MeV. We have included the finite widths 
of the D\ and T)|, and have set i] = 0. 
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(4.16) 


dr 

dcosct 


dr 


dcosa 


D* 


dr 

dcosct 


where 


dr 


dcosct 


dr 

dcosct 


3 sin 2 a 

(£„-A . 21 ) 2 + ri,./4’ 

An(l + 3 cos 2 a) - A s A d i2y^(l - 3 cos 2 a) cos r] sin 0 + 2M 2 sin 2 0 

(K-An^ + r^/d 


Here A 2 i and An are the resonant pion energies (averaged over charge states, A 2 i = 417 MeV 
and An = 383MeV), and 


A s - 

Adi — 1 + 


Vs f 


E-„ 


2 (fc/AJ El - ml ' 

8hi 


3 m c h 


(4.18) 


give the relative strength of the s-wave and d- wave transitions. Since A s contains a poten¬ 
tially large phase space enhancement, we keep A s sin 0 in our expressions. The two terms 


in Eq. ( 4.16 ) are normalized correctly with respect to each other, but the overall normal¬ 
ization is arbitrary. Note that while the resonances overlap, they do not interfere in this 
distribution, a feature which follows directly from Lorentz invariance. In Fig. [lj we show the 
distribution in cos a for E n = (380,400,420) MeV, for several scenarios of D\ decay: pure 
d- wave {A s = 0 and '0 = 0 ), pure s-wave {Adi = 0 and 0 = 7 t/ 2), and a mixed case with 
h\ = 0 and 0 = 9°, and |/| = 0.85. This latter case (actually, two cases, with / = ±0.85) 
corresponds to the situation in which there is no enhancement of the d-wave width, and the 
s-wave width is adjusted to give the correct total width of the D We have also chosen a 
value of 0 taken from the quark model [F|, and have set 77 = 0. We see that it should be 


possible to distinguish between these various scenarios. 


2. One-pion distributions 


We may also consider another distribution which is tied more closely to the fragmentation 
process by which the excited heavy meson is initially produced. In the decay {D^D^) —> 
H* 7 r, let 6 be the angle between the momentum of the pion and the fragmentation axis, as 
measured in the excited meson rest frame. In this frame, the fragmentation axis points back 
to the hard event in which the heavy quark was initially produced. The angular distribution 
in cos 6 depends not only on the quantities {h,hi,f, 77 , 0 }, but on the helicity distribution 
parameter W 3/2 [|HJ. This parameter describes the alignment with which the light degrees 
of freedom of Jn — \ are produced in the creation of the D\ or D\. The probabilities of the 
various helicity states along the fragmentation axis are given by 


nfi = D = nj? = -i)= ±3/2, 

P(J? = i) = P(Jf = -i)= |(1 - w 3/2 ). (4.19) 
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FIG. 2. The differential distribution in cos0, (a) without, and (b) with the finite widths of the 
D\ and D 2 . We show curves for E n = (380,400,420) MeV, identified as in Fig. 1, and have set 
= 9°, / = -0.85, hi = w 3 /2 = T) = 0. 


The parameter w 3 / 2 is a nonperturbative parameter of QCD, which is well defined only in 
the heavy quark limit. In Ref. [ T5| , data from ARGUS [H)| on the decay D 2 —■► Dn was 
used to set the 90% confidence level upper limit W 3/2 < 0.24. Models based on perturbative 
QCD have yielded an estimate w 3 / 2 ~ 0.25 l[l7|.p~8l . 

The distribution in cos 9 is considerably more complicated when hnite width effects are 
included, in part, because the interference between the D\ and D 2 resonances may not 
be neglected. (We note that when w 3 / 2 is extracted from the decay D* 2 —> Dn, there are 
no interference effects, since the decay Di —> Dn is prohibited.) After a straightforward 
calculation, we find 


dT 


dr 


oc 


d cos 9 d cos 9 


d T 


D* 


dcosd 


dr 


D 1 


dcosd 


(4.20) 


Di-D* 


where 


dr 


dcosd 

dr 


dcos# 


3(1 + cos 2 6 ) + 2u> 3 / 2 (1 — 3 cos 2 9) 

(£„-A 21 p + ri,./4 

1 


(4.21) 


[3(1 + cos 2 9) + 2w 3/2 (1 — 3 cos 2 9) 


Di (K-An ) 2 + r 2 Dl /4 

+ A s Adi^V2s\nijj cos7?(—1 + 2w 3 / 2 )(l — 3 cos 2 9) + 4 A 2 S sin 2 y;} 


dr 


dcos# 


(-1 + 2 ^ 3 / 2 )(1 - 3 cos 2 9) 


Di-D* 


2 m - A 21 )2 + r 2 /4] l(E n - Any + I%/4\ 


(4.22) 


x {{Adi ~ A s \[2 sin^cos rf) 4 (E n - A n )(£ ff - A 2i ) + T Dl V D * 

- A s 2\f2 sin^sin?y {E n T Dl - E n T D * + A U T D * - A 2i r Dl ) } . 

Once again, the three terms in Eq. (|4.2C|) are normalized only with respect to each other. 
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FIG. 3. The differential distribution in cos0, for W 3/2 = 0, in the case of (a) pure d -wave D\ 
decay, and (b) pure s-wave D\ decay. We have set / = —0.17, rj = 0, and E n = (380,400,420) MeV. 




FIG. 4. The differential distribution in cos 9, for W 3/2 = 0.25, in the case of (a) pure d- wave D\ 
decay, and (b) pure s-wave D\ decay. We have set / = —0.17, 77 = 0, and E n = (380,400,420) MeV. 

To explore the importance of including the finite widths of the excited resonances, in 
Fig. H we compare the distribution in cos 6 with (a) r(Z?i) = Y (D%) = 0 and (b) r(Z?i) = 
22 MeV, Y(D*) = 28 MeV, and scanning over E n = (380,400,420) MeV. For the purpose of 
illustration, we set W 3/2 = 0 and choose the parameters h\ = 0, ip = 9° and / = —0.85 for 
the D 1 decay. 

I 11 Fig. |3], we set w 3 / 2 = 0 and compare pure d -wave and pure s-wave D t decays, for 
/ = —0.17 and = (380,400,420) MeV. In Fig. fl] we do the same for w 3 / 2 = 0.25, in Fig. |] 
for W 3/2 = 0.75, and in Fig. |] for w 3 / 2 = 1.0. (For W 3/2 = 0.5, all the distributions are 
flat.) Together, Figs. give some idea of the sensitivity of the distribution in cos 9 to the 
various parameters describing the decay. Note that the sensitivity to W 3/2 is considerably 
enhanced if the D 1 decays via d- wave emission rather than s-wave emission. 
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FIG. 5. The differential distribution in cos 9, for W 3/2 = 0.75, in the case of (a) pure d -wave D\ 
decay, and (b) pure s-wave D\ decay. We have set / = —0.17, ?? = 0, and E n = (380,400,420) MeV. 




(a) (b) 

FIG. 6. The differential distribution in cos#, for rc 3 / 2 = 1-0, in the case of (a) pure d-wave D 1 
decay, and (b) pure s-wave D\ decay. We have set / = —0.17, i] = 0, and E n = (380,400,420) MeV. 

V. SUMMARY AND OUTLOOK 

Our detailed analysis of the leading effects of heavy spin symmetry violation on the 
properties of excited charmed and bottom mesons has led to a number of interesting results. 
In particular, the width of the D 1 , previously thought to be anomalously large, is seen 
actually to be of a natural size. Our predictions for the properties of excited bottom mesons 
agree well with the minimal data which exist so far, and will be tested soon in more detail. 
We presented detailed angular distributions for strong decays of excited charmed mesons, 
which will eventually provide a more stringent test of the predictions of our formalism. 

Excited heavy mesons are important both for their own sake and for the insight they 
give into the heavy quark expansion. Since the most accurate determinations of the CKM 
matrix element V c b involve the theoretical application of HQET, it is crucial to understand 
how well the m c , mb —■► 00 limit approximates the real world. It is a matter of more than 
academic interest whether the large width of the Di can be explained naturally within the 
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heavy quark expansion, because the answer to this question affects our willingness to trust 
that the charm quark may be treated as heavy in other contexts. Similarly, it is worthwhile 
to do one’s best to extract parameters such as w 3 / 2 , / and ^ from pion angular distributions 
in strong decays. Doing so, we learn not only about excited heavy mesons themselves, but 
about whether one can indeed explain their properties consistently in the context of the 
heavy quark expansion. 
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